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ABSTRACT 

We consider the backreaction of the magnetic field on the magnetic dynamo coefficients 
and the role of boundary conditions in interpreting whether numerical evidence for sup- 
pression is dynamical. If a uniform field in a periodic box serves as the initial condition 
for modeling the backreaction on the turbulent EMF, then the magnitude of the turbulent 
EMF and thus the dynamo coefficient a, have a stringent upper limit that depends on the 
magnetic Reynolds number Rm to a power of order —1. This is not a dynamic suppression 
but results just because of the imposed boundary conditions. In contrast, when mean field 
gradients are allowed within the simulation region, or non-periodic boundary are used, the 
upper limit is independent of Rm and takes its kinematic value. Thus only for simulations of 
the latter types could a measured suppression be the result of a dynamic backreaction. This 
is fundamental for understanding a long-standing controversy surrounding a suppression. 
Numerical simulations which do not allow any field gradients and invoke periodic boundary 
conditions appear to show a strong a suppression (e.g. Cattaneo & Hughes 1996). Sim- 
ulations of accretion discs which allow field gradients and allow free boundary conditions 
(Brandenburg & Donner 1997) suggest a dynamo a which is not suppressed by a power of 
Rm- Our results are consistent with both types of simulations. 

Subject Headings: magnetic fields; galaxies: magnetic fields; Sun: magnetic fields; 
stars: magnetic fields; turbulence; accretion discs. 
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1. Introduction 

A leading candidate to explain the origin of large scale magnetic fields in stars and 
galaxies is mean-field turbulent magnetic dynamo theory (Moffatt 1978; Parker 1979; Krause 
& Radler 1980; Zeldovich et al. 1983, Ruzmaikin et al. 1988, Beck et al. 1996). The theory 
appeals to a combination of helical turbulence (leading to the a effect), differential rotation 
(the Q effect), and turbulent diffusion to exponentiate an initial seed mean magnetic field. 
The total magnetic field is split into a mean component and a fluctuating component, and 
the rate of growth of the mean field is sought. The mean field grows on a length scale much 
larger than the outer scale of the turbulent velocity, with a growth time much larger than 
the eddy turnover time at the outer scale. A combination of kinetic and current helicity 
provides a statistical correlation of small scale loops favorable to exponential field growth. 
Turbulent diffusion is needed to redistribute the amplified mean field rapidly to ensure a net 
mean flux gain inside the system of interest. Rapid growth of the fluctuating field necessarily 
accompanies the mean-field dynamo. Its impact upon the growth of the mean field, and the 
impact of the mean field itself on its own growth are controversial. 

The controversy results because Lorentz forces from the growing magnetic field react 
back on and complicate the turbulent motions driving the field growth (e.g. Cowling 1959, 
Piddington 1981, Kulsrud & Anderson 1992). It is tricky to disentangle the back reaction of 
the mean field from that of the fluctuating field. Analytic studies and numerical simulations 
seem to disagree as to the extent to which the dynamo coefficients are suppressed by the 
back reaction of the mean field. Some numerical studies (e.g. |Cattaneo fc Vainshtein ( 1991 )| , 



Vainshtein fc Cattaneo ( 1992 )| , |Cattaneo ( 1994 )| , |Cattaneo fc Hughes (1996]|) and analytic 



studies (|Gruzinov fc Diamond ( 1994 )| , |Bhattacharjee fc Yuan ( 1995 )| , [Kleeorin et al., (1995)| ) 



argue that the suppression of a takes the form a ~ cr°)/(l + R P M B 2 /vq) where a^ ^ is the 
value of a in the absence of a mean field, Rm is the magnetic Reynolds number, B is the 
mean field in velocity units, t>o is the rms turbulent velocity, and p is a number of order 1. 
Such a strong dependence on Rm would prevent astrophysical dynamos from working, as 
Rm is usually ^> 1. 

Other numerical studies (Brandenburg & Donner 1997) and analytic studies (e.g. Kraich- 
nan 1979; Field et al, 1999, Chou & Fish 1999) suggest that p = 0, so a ~ at® /{I + B 2 /v%) 
in the fully dynamic regime. In particular Field et al. (1999) considered an expansion in the 
mean magnetic field (see also Vainshtein & Kitchatinov 1983; Montgomery & Chen 1984; 
Blackman & Chou 1997), and were able to derive the effect of the nonlinear back reaction 
on a in the case for which VB = 0. Their result is expressed in terms only of the difference 
between the zeroth-order kinetic and current helicities. They find that Rm does not enter 
strongly, except possibly by suppressing the difference between the zeroth-order helicities, 
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an effect which cannot depend upon B and is not therein constrained. 

Blackman & Field (1999) have shown that some of the analytic approaches (e.g. Bhat- 
tacharjee & Yuan 1995; Gruzinov & Diamond 1994) which employ a use of Ohm's law dotted 
with the fluctuating component of the magnetic field, do not distinguish between turbulent 
quantities of the base (zeroth-order) state and quantities which are of higher order in the 
mean field. This distinction is important. When it is made, many arguments for suppression 
via such approaches fall through. Note that Blackman & Field (1999) do not prove that the 
dynamo survives back reaction as a result of their considerations, only that some analytic 
approaches to the problem can be challenged. 

Despite this challenge, the apparent result of extreme a suppression is seen in the simu- 
lation of Cattaneo & Hughes (1996). These authors externally force the turbulence, imposing 
periodic boundary conditions and a uniform mean field, and find that the suppression of a 
involves Rm in the form given above. By contrast, the simulation of Brandenburg & Donner 
(1997) suggests that ana-0 dynamo may in fact be operative in an accretion disc whose 
turbulence is self-generated by a shearing instability, without Rm entering the suppression. 
The latter simulation does not employ periodic boundary conditions and allows gradients in 
mean fields. 

In this paper we show that the suppression of a depends crucially on the boundary 
conditions. We find that the mean quantities are defined by averaging over a periodic box a 
has an upper limit that depends on a factor of R~^ , with p ~ 1. In the presence of mean field 
gradients and non-periodic boundary conditions, however, we find that the upper limit on 
the dynamo coefficients is significantly larger, and Rm is not involved. The small upper limit 
in the periodic box case does not represent a dynamical suppression but rather an apparent 
suppression which is just a results from the boundary conditions. The results herein may be 
a step toward resolving controversies surrounding numerical suppression experiments. 

Central to the discussion is the equation for the time evolution of magnetic helicity. This 
equation was also employed by Seehafer (1994), who derived a suppression of a apparently 
consistent with that of Keinigs (1983) (and qualitatively consistent with the Cattaneo & 
Hughes (1996) simulation). The techniques of these additional two papers are different and 
one should note that they do not separate zeroth from higher order quantities. 

Section 2 reviews the basic formalism of the dynano coefficient expansion in orders of 
B. Section 3 shows that constraints on the magnitude of the EMF (and thus a) results 
from dotting Ohm's law for the fluctuating electric field with the fluctuating magnetic field, 
taking the average and expanding to second order in the mean magnetic field. The results 
depend on the boundary conditions. For a periodic box, the upper limit on a is too small 
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for a dynamo to work in practice, but this does not represent a dynamical suppression. In 
section 4 we interpret the results in terms of helicity flow and we discuss implications with 
respect to previous studies. Section 5 is the conclusion. 

2. Basic Formalism 

The basic formalism employed herein is discussed in Field et al. (1999) and Blackman 
& Field (1999). The formalism combines some aspects of the standard textbook treatment 
(e.g. Moffatt 1978) with the modification that fluctuating quantities are divided into (1) 
a zeroth-order turbulent base state whose correlations are homogeneous, stationary and 
isotropic (though not necessarily mirror-symmetric!), and (2) a contribution which depends 
on the presence of a non-zero mean field (see also Vainshtein & Kitchatinov 1983; Kitchatinov 
et al. 1994). This higher-order contribution is definitely not isotropic and not necessarily 
homogeneous or stationary. 

More specifically, the induction equation describing the magnetic field evolution is 

d t B = V x (V x B) + AV 2 B, (1) 

where A = r]c 2 /Ait is the magnetic viscosity, corresponding to resistivity rj. Here B = b + B 
is the magnetic field in velocity units, obtained by dividing by V^rp, and b and B are the 
fluctuating and mean components of B respectively. We assume incompressibility. 

The equation for the mean field derived by averaging ([]]) is 

d t B = V x (v x b) - V • VB + B ■ VV + AV 2 B. (2) 
The equation for b is given by subtracting (§) from ([!]), which gives 

d t b = V x (v x b) — V x (v x b) + V x (v x B) + V x (V x b). (3) 

The term B ■ VV in (|2|) describes the f2-effect of differential rotation, and will not be 
discussed further here, while the term V • VB can be eliminated by changing the frame of 
reference to one moving with V; both terms will be ignored in what follows. The dynamo 
theorist must find the dependence of the turbulent EMF (v x b) on B so that (|2]) can be 
solved. 

In the absence of the mean velocity fields, all mean vectors can be written in terms of 
the mean magnetic field. In particular, we have 

(v x b) = OijBj - P ijk djB k + j m O(B/R 2 ) + (4) 

where ct^, (3^^ and 7^ are explicit functions of correlations of turbulent quantities, but can 
depend implicitly on B (Moffatt 1978) through their dependence on the induction equation 
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for the fluctuating field. The order at which there is no implicit dependence on B is the 
zeroth-order base state (see Field et al 1999). The expansion order parameter is | B | / 1 1>^ ^ | , 
which is indeed << 1 for the early dynamo evolution, and < 1 in the Galaxy at present. 
In particular, we have b = b^ ^ + J2n^ n \ and similarly for v, where S n b( n ) < b^) and 
n indicates the number of powers of |B|/|b|. The zeroth order base state correlations are 
composed of products of and v^ ) and have no dependence on the mean field. The zeroth 
order base state is taken to be homogeneous and isotropic-the violation of isotropy comes 
from the contributions due to higher order fluctuating quantities, whose isotropy is broken 
by the mean field. Note that the zeroth order state need not be reflection invariant, and it 
is important for dynamo theory that it is not. 

Correlations between higher order quantities can be reduced to correlations of zeroth 
order quantities times the respective products of n linear functions of B. Thus for example, 
b^ 2 ) is the anisotropic component of the fluctuating magntetic field which depends on two 
powers of B, and is found by twice iterating terms like b ■ B in the induction equation to 
obtain an approximate solution in terms of b^ -* and v^ -*. 

To zeroth order, the a tensor can be written 

al? = « ( %, (5) 

which highlights the isotropy of this zeroth-order quantity. In our previous work (Field 
et al. 1999) we have used the induction equation for the fluctuating components of the 
magnetic field and the Navier-Stokes equation for the fluctuating velocity to find the form 
of a in terms of correlations of the zeroth-order products (see also Blackman & Chou 1997). 
Calculating the turbulent EMF in the absence of gradients of B, to first order in B, gives 
(v x b)( x ) = a^B, where is the sum of kinetic and current helicities associated with the 
zeroth order state, namely 



« (0) = 4 



v(°)(t)/-V x ^(t')dt') - (b(°)W/-V x b(°)(t')df 



_|t c f( v (o) . V x v<°>) - (b(°) • V x b(°) 



(6) 



where v is the turbulent velocity and t c defined as the correlation time of the scale of 
the turbulence which dominates the averaged quantity. If we adopt a Kolmogorov energy 
spectrum (i.e. kb%., hv\ oc k 2 Ek oc A; 1 / 3 ), then it might appear that the dominant contributions 
to the terms of @ come from large k. However, Pouquet et al. (1976) showed that if the 
forcing is at the outer wavenumber k = most of the energy and helicity is concentrated 
there, and the turbulence for k > 3ko is locked up into Alfven waves which do not contribute 
to correlations. It is therefore likely reasonable to assume that any helicity in the zeroth-order 
state is concentrated near ko, in which case t c ~ (fo^o) -1 ~ L/vq. 
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The first term in (|]) was first derived by Steenbeck, Krause, & Radler (1966). The 
second, current helicity, term in (|6|) was first derived by Pouquet et al. (1976); neither paper 
made the necessary distinction between zeroth and higher-order quantities. 

In the next sections we will not re-derive the form of a/°) in terms of b^ ^ and v^ ); instead 
we will derive an independent upper limit on from the use of Ohm's law, the definition 
of the electric field in terms of the vector potential, and the equation for magnetic helicity 
evolution. We will invoke the assumption that the zeroth-order base state is isotropic and 
homogeneous, and we will assume that that all anisotropies and inhomogeneities of higher- 
order correlations are due to mean fields. 



We will need the Reynolds relations QRadler (1980)|) , i.e. that derivatives with respect 



to x or t obey dt, x (XiXj) = (dt, x (XiXj)) and (XiXj) = where X; = Xi + Xi are components 
of vector functions of x and t. For statistical ensemble means, these hold when correlation 
times are small compared to the times over which mean quantities vary. For spatial means, 
defined by pQ(x, t)) = V' 1 j -X»(x + s, t)ds, the relations hold when the average is over a 
large enough V that L <C Vs < R < D, where D is the size of the system R is the scale 
of mean field variation and L is the outer scale of the turbulence. Note that the scale of 
averaging is less than the overall system size. 

Constraints on the tubulent EMF for periodic and non-periodic boundary conditions 

a. Constraint equations 

Let the electric field E, like B, be divided into a mean component E and a fluctuating 
component e. Ohm's law for the mean field is thus 

E= (-c _1 V x B + V j) = -c -1 (v x b) + r/ J (7) 

for the case V = 0, where J is the current density and rj is the resistivity. We also have 

(E ■ B) = E ■ B + (e ■ b) = -cT 1 (v x b) ■ B + r]3 ■ B + (e • b) (8) 

where we have used (0). 

A second expression for (E-B) also follows from Ohm's law without first splitting into 
mean and fluctuating components, that is 

(E-B) = (-c 4 (VxB).B + i)J.Bj 

= r? ( J • B) = rj J • B + 7] (j • b) = 7]3 ■ B 

+c~ 1 X (b • V x b) . (9) 

By substituting (||) into (|[), we obtain 

(T 1 (v x b) • B = -c _1 A (b ■ V x b) + (e ■ b), (10) 
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an equation which will now constrain (v x b). However, we must expand ([R]) to second 
order in B (as defined in section 2) to constrain the turbulent EMF (v x b). This is because 
to zeroth order, the left hand side of ( p!U|) vanishes directly. To first order, the left side would 
be (v x b)^B, but (v x b)(°) = 0, since vector averages of zeroth order quantities vanish. 
To second order in B then, fllTf ) implies that 

c- 1 (v x b) (1) • B = -c _1 A (b • V x b) (2) + (e ■ b) (2) . (11) 

Because Rm >> 1, significant limits on (v x h)^ and thus on come from the (e • h)^ 
term above. The result of Seehafer (1994) and Keinigs (1983) amount to the ([11]) with the 
last term equal zero, but without distinguishing the order in mean fields (i.e. without the 
superscripts). We now focus on this last, term keeping in mind that it is second order in 
mean fields. 

Since (e • b)( 2 ) is second order in B, its most general form will be expressible as a sum 
of terms which each involve products of two types of quantities: 1. correlations of scalar 
or pseudoscalar products of zeroth order quantities and 2. quadratic scalar or pseudoscalar 
functions of B. Now (e • b) can be written as a sum a time derivative and spatial divergence. 
Consider e in terms of the vector and scalar potentials a and 4>: 

e = -V0 - (l/c)<9 t a. (12) 

Dotting with b = V x a and averaging we have 

(e ■ b) = -(V0 ■ b) - (l/c)(b ■ <9 t a>. (13) 

After straightforward algebraic manipulation, application of Reynolds rules and V • b = 0, 
this equation implies 

(e ■ b) = -(1/2)V ■ (0b) + (1/2)V -(axe)- (l/2c)Sfc(a ■ b) = -d h° + dji = d^V, (14) 

where we have defined a helicity density 4-vector for fluctuating quantities 

[ho, hi) = [(l/2c)a • b, (1/2) ■ (0b) - (l/2)(a x e)], (15) 

and the overbar is used, as always, to mean the same thing as the brackets. 

b. Constraints for periodic boundary conditions 

We now investigate the implications of ( |14"D for simulations of type performed by Catta- 
neo & Hughes (1996), where the brackets brackets are interpreted as a spatial average over 
a periodic box. Under these conditions, there are two important consequences. First, note 
that the second two terms of (|i~4|) vanish upon conversion to surface integrals and we have 



(e-b) = -(l/2c)0i(a-b), 



(16) 
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which is gauge invariant. The second consequence of the periodic box is that d t B = for 
incompressible flows. This follows simply from (0): the last three terms of (0) would vanish 
as they are all surface integrals. Using Reynolds rules and vector identities, the second 
term can be written [V x (v x b)]j = (di(biVj)) — (di(vibj)) which also vanishes by surface 
integration. 

The two consequences just discussed can be used to show that fll6f ) vanishes for a periodic 
box, and thus the only contribution to the right of (|TT|) will come from the first term on the 
right. To second order in mean quantities, assuming b(t = 0) = and that all times are far 
enough from t = such that b(t) does not correlate with any finite a(0), we have 

(a-b)^ = J (d t ,a.W(t')-b(°\t))dt'+ J (dt^it'yj d t »b^(t"))dt'dt"+ J (a^(t)-d t/ b^(t')}dt'. 

(17) 

To express (|T7|) explicitly in terms of mean fields, we use the equations of motion for b 
and a. The use of d t b from (§) for the last two terms of ([[7]) leads directly to contributions 
depending on products of the mean fields B or V and turbulent quantities b and v. Consider 
now the equation for a which comes from uncurling the equation for b, namely 

d t SL = (v x b) - (v x b) + (v x B) + (V x b) + V0, (18) 

where 9 is an arbitrary scalar field. When ( |T8"D is used in QT7| ) in the first and second terms on 
the right of ([T7|), the periodic box nullifies the contribution from W. All other contributions 
depend only on products of v, b, B and V. Thus when V = 0, the only remaining mean 
field is B. Thus for a periodic box, (a • b)( 2 ) must be second order in B. Then, when 



plugged into ( [Lq) the time derivative will act on some quadratic function of B multiplied 
by correlations of zeroth order. Since the zeroth order quantities are time independent, 
isotropic, and homogeneous, the function of B must be a scalar, denoted F, and we have 

d t (a ■ b)& = d t (F(B)^Q^) = gi 0) ^(F(B)( 2 )) = 0, (19) 
where Q± is a scalar or pseudoscalar correlation of zeroth order quantities. The last equality 



of (p"9D follows form stationarity of zeroth order quantities, and our proof that B is time 
independent over the time scales of interest for the periodic box. We therefore conclude that 
<9t(a • b) = in fllEp . This result relates to the the fact that for a periodic box, there is no 
periodic mean vector field A whose curl is everywhere equal to B. The divergence of B is 
still equal to zero, so Maxwell's equations are satisfied, but B is the only non-trivial mean 
field. 

Since in the Cattaneo & Hughes (1996) simulation B =constant in both space and time, 
(v x b)( 2 ) = a^B 2 . Using this, and (|1|), (0) and ([II]), we obtain 

= 

B 
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Field et al. (1999) showed that conclusions about cr°) are also conclusions about a to all 
orders, by relating the fully non-linear a to and showing that in the limit of large B, a 
is not catastrophically affected when B is large. Thus is an upper limit to a, and so the 
result ( |2"0D shows that a will be small when the brackets indicate an average over a periodic 
box. The important point is that this is not a dynamical suppression from the backreaction 
but a constraint on the magnitude of which is imposed by the boundary conditions. 
Notice that it is a constraint on the zeroth order quantity, and so it cannot represent the 
effect of backreaction. 

c. Constraints for non-periodic boundary conditions 

If the averaging brackets are not over a periodic box, or if the scale of the averaging 
is << than the overall scale size of the system, then the divergence terms in (14) do not 



vanish. In addition, the V0 term in (|i~8|) will contribute to ([161) . In this case, each term 



on the right of (14) is not gauge invariant. Thus, the only constraint we can make on the 



magnitude of the right side of (II) is on the sum of all the terms together. Writing down all 



possible second order terms up to one spatial derivative in B, we have 

c(e • b>^ = Qf B 2 + gf B • V x B + 0(B~ 2 L 2 /R 2 ) + (21) 

where and are correlations of zeroth order averages, L is the outer turbulent scale 
and R is the mean field variation scale. The quantity Q2 must have units of velocity, and 
thus be of maximum order vq since it depends only on turbulent quantities. The quantity 
must have dimensions of viscosity, and must of maximum ~ vqL since it too depends 



only on turbulent quantities. The combination of terms in (|21|) is the same form of the 
combination of terms entering on the left side of ( |10D which would result from using (|4]). 
That is, since (v x b) • B = a^B 2 - (3^B • V x B, we have = a<® and Qf = -/3<®. 
Thus for simulations in which the mean values are not taken over a periodic box, there is 
certainly no a priori restriction on a^°\ Since now a < a^°\ any simulation result indicating 
suppression on a under these relaxed boundary conditions would indeed be a dynamical 
suppression. So far there are no simulations which invoke such boundary conditions that 
show catastrophic suppression (c.f. Brandenburg & Donner 1997). 

4. Discussion 

Section 3 shows that periodic boundary conditions impose an upper limit on a that 
does not represent a dynamical suppression. Non-periodic boundary conditions or a finite 
scale separation between system size and mean field gradient scale allow for a much higher 
upper limit on a, namely the kinematic limit. The dynamical backreaction is testable only 
in simulations of the latter type. 

4.1 Relation to magnetic helicity 
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Here we point out a connection to magnetic helicity. Repeating Eqns. (0), ([Tq ) and 
(0) for the total E and B gives 

(E • B) = E • B + (e • b) = ±dpH» = Uji^ + -fyF ~ 0, (22) 

where if M is the total magnetic helicity 4-vector (Field 1986) defined exactly as in QT5D but 
with all fluctuating quantities replaced by their total values. Similarly, is the helicity 4- 
vector associated with the mean fields. The last similarity in ( |22"D follows because Ru >> 1 
(A ~ 0) in the astrophysical plasmas of interest. Using (e • b) = d^h , Eqn. (p2|) then shows 
that any non-negligible (e ■ b) requires a finite but opposite d^h^. 

In general, for non-vanishing turbulent EMF, the (e • b) must be non-zero, and thus 
the 4-divergences in fl2"2]) cannot vanish separately. Interestingly, when (|2"2"D is integrated 



over the total volume inside and outside of the rotator, and interpreted in terms of a flow 
of relative magnetic helicity (Berger & Field 1984), it can be shown that a working dynamo 
implies an associated magnetic energy flow through the magnetic rotator of interest which 
likely leads to an active corona (Field & Blackman 2000; Blackman & Field 2000). 

4.2 Implications for Previous and Future Studies 

The use of periodic boundary conditions in simulations appears to be unsuitable for 
testing the suppression of a in a real dynamo unless the scale of mean field variations is 
much smaller than the scale of the periodicity. If periodic boundary conditions are used, 
one must also be careful about causality issues. The scale separation should at minimum 
be large enough such that the Alfven crossing time across the box is longer than correlation 
time of the fluctuating quantities, and possibly even longer than the time scale for mean 
field variation. Thus, the box could be periodic, but the dynamics of interest would occur 
in a non-periodic sub-region. The brackets which we have used to indicate averages, would 
thus represent averages over this sub-region, not the entire volume. Alternatively, the box 
could be non-periodic. 

The numerical simulations of Cattaneo & Hughes (1996) do not allow for any mean 
field gradients and employ periodic boundary conditions. The strong a reduction seen there 
is consistent with our suggestion that the suppression may not be dynamical, but may 
instead be a result of the boundary conditions. In contrast, the shearing box accretion disk 
simulations of Brandenburg & Donner (1997) do employ non-periodic boundary conditions 
and allow mean field gradients. Interestingly, they do find that something like a mean-field 
dynamo is operating therein. The limited suppression that they find does not involve Rm- 

5. Conclusions 



We have suggested that the cause for apparent a suppression in numerical simulations 
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which use periodic boundary conditions may not result from dynamics, but from rather from 
a choice of boundary conditions. If the boundary conditions enforce all mean field gradients 
and spatial divergences to vanish, then the upper limit on a is given by fl2H|). For non- 
periodic boundary conditions or a box with significant scale separation between the mean 
field and box size, the upper limit on the turbulent EMF is given by the kinematic value. 
This would be a consistent interpretation of the large suppression reported by Cattaneo 
& Hughes (1996). In contrast, Brandenburg and Donner (1997) interpret disk simulations 
which use non-periodic boundary conditions and do not find such a strong suppression. In 
summary, our results are consistent with seemingly contradictory simulations. 

Working dynamos in real astrophysical bodies (even in the kinematic approximation) 
require mean field gradients and scale separations between the overall system scale, mean 
field averaging scale, and fluctuating scale. In order to disentangle boundary effects from 
dynamical ones, future simulations of a suppression should include non-periodic boundary 
conditions or allow the mean field to change over scales smaller than the size of the overall 
box. This is a challenging task. 
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REFERENCES 

Beck, R. et al. 1996, Ann. Rev. Astron. Astrophys., 34, 155 
Berger, M. A. & Field, G. B., 1984, J. Fluid Mech. 147, 133. 
Bhattacharjee, A. & Yuan Y., 1995, ApJ, 449, 739. 
Blackman, E.G. & Chou, T., 1997, ApJ 489 L95 
Blackman, E.G. & Field, G.B., 1999, ApJ, in press. 
Blackman, E.G. & Field, G.B., 2000, MNRAS, submitted. 
Brandenburg, A. & Donner, K.J., 1997, MNRAS 288, L29. 
Cattaneo, F. 1994, ApJ, 434, 200 

Cattaneo, F. & D.W. Hughes, 1996, Phys. Rev. E., 54, 4532. 
Cattaneo, F. & Vainshtein, S.I., 1991, ApJ, 376, L21. 



- 12 - 



Chou, H.S. & Fish, V., 1999, ApJ, submitted. 

Cowling, T.G., 1957, Magnetohydrodynamics, (New York: Interscience) . 

Field, G. B., 1986, Magnetic Helicity in Astrophysics, in Magneto spheric Phenomena in 
Astrophysics, R. I. Epstein & W. C. Feldman, eds., AIP Conference Proceedings 144 
(Los Alamos: Los Alamos Scientific Laboratory, 1986), 324. 

Field, G. B., 1995, in Proceedings of the International Conference on Plasma Physics, ed.P. 
H. Sakanaka & E. del Bosco (Woodbury, New York: AIP Press), 416. 

Field G.B. & Blackman E.G. 1999, in Highly Energetic Physical Processes and Mech- 
anism for Emis son from Astrophysical Plasmas, IAU Symp. 195 eds. P.C.H. Martens 
and S. Tsuruta. 

Field, G. B. & Blackman, E. G. & Chou, H., 1999, ApJ, 513 638. 
Gruzinov, A. & Diamond PH., 1994, PRL, 72, 1651. 

Kitchatinov, L. L., Pipin V.V. and Riidiger, G., & Kuker, M. 1994, Astron. Nachr., 315, 
157 

Keinigs, R.K., 1983, Phys. Flu. 26 2558. 

Kleeorin, N., Rogachevskii, & Ruzmaikin, A., 1995, A&A 297 159 
Kleeorin, N. & Ruzmaikin, A., 1982, Magnetohydrodynamics, N2, 17. 
Kraichnan, R. H. 1979, Phys. Rev., 113, 1181 
Kraichnan, R.H, 1979, PRL, 42, 1667 

Krause, F., & Radler, K.-H. 1980, Mean-Field Magnetohydrodynamics and Dynamo Theory, 
Oxford: Pergamon Press 

Kulsrud, R. M., & Anderson, S. W. 1992, ApJ, 396, 606 

Moffatt, H. K. 1978, Magnetic Field Generation in Electrically Conducting Fluids, Cam- 
bridge: Cambridge University Press 

Montgomery, D. & Chen, H., 1984, Plasma Physics & Controlled Fusion, 26, 1189. 

Parker, E. N. 1979, Cosmical Magnetic Fields, Oxford: Clarendon Press 

Piddington, J. H. 1981, Cosmical Electrodynamics, Malbar: Krieger 



-13- 



Pouquet, A., Frisch, U., and Leorat, J. 1976, JFM 77, 321 
Radler, K.-H. 1980, Astron. Nachr., 301, 101 

Ruzmaikin, A. A., Shukurov, A. M., and Sokoloff, D. D. 1988, Magnetic Fields of Galaxies, 
Dordrecht: Kluver Press 

Seehafer N., 1994, Europhys. Lett., 27 353. 

Steenbeck, M., Krause, F., & Radler, K.H., 1966, Z. Naturforsch., 21a, 369. 

Vainshtein S., 1998, PRL 80, 4879. 

Vainshtein, S.I. k Cattaneo F., 1992, ApJ, 393, 165. 

Vainshtein, S. I., and Kitchatinov, L. L. 1983, Geophys. Ap. Fluid. Dyn., 24, 273 

Zeldovich, YA. B., Ruzmaikin, A. A., and Sokoloff, D. D. 1983, Magnetic Fields in Astro- 
physics , New York: Gordon and Breach 



This preprint was prepared with the AAS IATeX macros v4.0. 



